Observation of CP Violation in Kl — > 7r + 7r e + e Decays 
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We report the first observation of a manifestly CP violating effect in the Kl — > 7r + 7r _ e + e _ decay mode. A large 
asymmetry was observed in the distribution of these decays in the CP-odd and T-odd angle (f> between the decay 
planes of the e + e _ and 7T + 7r~ pairs in the Kt, center of mass system. After acceptance corrections, the overall asym- 
metry is found to be 13.6 ± 2.5 (stat) ± 1.2 (syst)%. This is the largest CP-violating effect yet observed integrating 
over the entire phase space of a mode and the first such effect observed in an angular variable. 



PACS numbers: 13.20.Eb, 13.25.Es, 13.40.Ag, 14.40.Ag 

The KTeV E799 experiment at Fermi National Ac- 
celerator Laboratory recently reported the first observa- 
tion PJ of the four body decay mode Kt, — » tt + tt~ e + e~ . 
Based on 2% of the data, a branching ratio of 3.2 ± 
0.6 (stat) ± 0.4 (syst) x 10~ 7 was measured. In this 
paper, we report an analysis of the entire KTeV E799 
data from which the Kl — > 7r + 7r~e + e~ signal (shown in 
Fig. [I]) of 1811 events above background was obtained 
after the analysis cuts described below. We observed in 
these Kl — > ir + ii~e + e~ data a CP-violating asymmetry 
in the CP- and T-odd variable sin 6 cos 6, 
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in the Kl cms. This asymmetry implies, with mild as- 
sumptions, time reversal symmetry violation as well. The 
quantity sin^cosi^ is given by (h ee x n^n) • z(n ee ■ h^n), 
where the n' s are the unit normals and z is the unit vec- 
tor in the direction of the tttt in the Kl cms. 



FIG. 1. M^+j-e+j- invariant mass for events passing cuts. 

The observed asymmetry sin <f> cos <j> shown in Fig. |^ 
was 23.3 ± 2.3 (stat)% before corrections. Inspection of 
Fig. ^| shows that the asymmetry between the bins near 
sin (j> cos <fi — ±0.5 is considerably larger. As discussed 
below, this cannot be explained by asymmetries due to 
either the spectrometer acceptance or detector elements. 
Using the model of Ref. to correct for regions of 

Kl —> 7r + ?T~e + e~ phase space outside the acceptance of 



1 



the KTeV spectrometer (which have small asymmetry), 
an asymmetry integrated over the entire phase space of 
the Kl — > 7r + 7r~e + e~ mode of 13.6 ± 2.5 (stat)% was 
obtained, the largest such CP- (and T-) violating effect 
yet observed. In comparison, CPLEAR recently reported 
a 0.66 ±0.13 (stat)% T-violating asymmetry || between 

K° — » K and K — > K° transition rates. 
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FIG. 2. a) Observed and b) sin cos angular distribu- 
tions: The data are shown as dots. The histogram is a Monte 
Carlo simulation based on the model of Ref. [3]. 

The Kl — » 7r + 7r~e + e~ data were accumulated during 
the ten weeks of E799 operation. A proton beam with 
intensity typically in the range 3.0 - 3.5 x 10 12 protons 
per 23 second spill incident at an angle of 4.8 mr on a 
BeO target was employed to produce two nearly paral- 
lel K L beams for E799. The configuration of the KTeV 
E799 spectrometer consisted of a vacuum decay region, a 
magnetic spectrometer with four drift chambers, photon 
vetoes, eight transition radiation chambers, a Csl elec- 
tromagnetic calorimeter, and a muon detector. A total 
of 2.7 x 10 11 Kl decays were accumlated during the E799 
run. Details of the KTeV detector are given in Ref. 0]. 

The KTeV four track trigger §J selected 1.3 x 10 8 
events. Candidate Kl — » 7r + 7r _ e + e _ events were ex- 
tracted from these triggers by requiring events with four 
tracks that passed track quality cuts and had a com- 
mon vertex with a good vertex x 2 . To be designated 
as e 1 * 1 , two of the tracks were required to have opposite 
charges and 0.95 < E/p < 1.05 where E was the energy 
deposited by the track in the Csl and p was the momen- 
tum obtained from magnetic deflection. To be consistent 
with a pair, the other two tracks were required to 
have E/p < 0.90 and opposite charges. 

To reduce backgrounds arising from other types of Kl 
decays in which decay products have been missed, the 
candidate 7r + 7r~e + e~ were required to have transverse 
momentum P t 2 of the four tracks relative to the direction 
of the K L be less than 0.6 x 10" 4 GeV 2 /c 2 . This cut was 
91.8% efficient for retaining Kl — * 7r + 7r~e + e _ . 

The major background to the Kl — > 7r + 7r~e + e~ mode 
was Kl — > 7r + 7r~7r^, where -k° d was a Dalitz decay, 
7t° — > 7e + e~, in which the photon was not observed in 
the Csl calorimeter or the photon vetos. To reduce this 
background, all Kl — * 7r + 7r~e + e~ candidate events were 



interpreted as Kl — > tt + tt^tt q d decays. Under this as- 
sumption, the momentum squared P 2 of the assumed 
7T can be calculated in the frame in which the momen- 
tum of 7r + 7r _ is transverse to the Kl direction. P 2 was 



mostly greater than zero for Kl 
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decays ex- 



cept for cases where finite detector resolution produces a 
P%o < 0. In contrast, most of the Kl — > tt + it~ e + e~ de- 
cays had P 2 < 0. The requirement that all 7r + 7r~e + e~ 
had (P w o) 2 < -0.00625 GeV 2 /c 2 minimized the K L -» 
7r^7r — 7r^j background while retaining 94.8% of the signal. 

Other backgrounds were relatively minor. The largest 
of these was due to Kl — > Tr + n~j decays in which 
the photon converted in the material of the spectrom- 
eter. These events, which reconstructed to the Kl 



mass and survived the P t and P 2 cuts, were elimi- 
nated by requiring M e+e - > 2.0 MeV/c 2 . The M e+e - 
cut retained 95.3% of the Kl — > 7r + 7r~e + e - events. A 
third background due to accidental coincidence of two 
Kl — > 7r ± e T ^ decays (K e3 ) whose decay vertices over- 
lap was minimized by track and vertex \ 2 cuts. A fourth 
background due to S° — > Kit'jj where the proton from the 
A decay was misidentified as a ir + was made negligible by 
K° momentum and vertex \ 2 cuts. Finally, a fifth back- 
ground due to Ks — > 7r + 7r~e + e - decays was eliminated 
by requiring the energy of the irnee be <200 GeV/c 2 . 

The final requirement of the Kl — > 7r + 7r~e + e~ events 
was 492 MeV/c 2 < M„ ee < 504 MeV/c 2 . The magni- 
tude of the background under the Kl peak was deter- 
mined by a fit to the 7r + 7r _ e + e _ mass distribution out- 
side the signal region. From this fit, a Kl — * Tr + n~e + e~ 
signal of 1811 ± 43 (stat) events above a background of 
45 ± 11 events was obtained in the signal region. The 
45 event background was composed of residual Kl — > 
7r + 7r~7T£ (36 events), Kl — > tt + tt~j (4.0 events), over- 
lapping K e3 (3.5 events), cascade decays (1.3 events) and 
Kg — > 7r + 7r~e + e~ (0.2 events). 

Possible sources of false asymmetries were considered, 
including those due to backgrounds and asymmetries in 
the detector. To check for detector asymmetries, the co- 
pious Kl —* Tt^ir^it^ decay mode, which has a similar 
topology to Kl — ► 7r + 7r~e + e - except for the presence of 
an extra photon in the Csl, was used. This mode is ex- 
pected to have no asymmetry in the (f> distribution formed 



using the n + n 



In a sample of approximately five 



million Dalitz decays, an asymmetry of —0.02 ± 0.05% 
was observed. The small background under the Kl was 
determined not to contribute significantly to the asym- 
metry in the Kl mass region since the asymmetry of 
the sideband regions below and above the Kl mass was 
measured to be 3.1 ± 5.1% and —2.3 ± 9.2% respectively. 

To perform an acceptance correction for loss of events 
due to spectrometer geometry, trigger, reconstruction 
efficiency, and analysis cuts, we modeled the Kl — > 
7r + 7r~e + e~ decays and simulated the response of the 
KTeV detector elements. The Kl — * 7r + 7r~e + e~ decay 
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mode is expected [|2|-|4| to proceed via both CP-violating 
and conserving amplitudes and exhibit both direct and 
indirect CP-violation. The dominant CP-violating am- 
plitude is indirect and proceeds via an initial decay of 
the Kl into 7r + 7r~ followed by one of the pions undergo- 
ing inner bremsstrahlung with the resulting photon in- 
ternally converting to an e + e~ pair. The dominant CP- 
conserving amplitude is the emission of an Ml photon 
at the 7r + 7r~ decay vertex followed by internal conver- 
sion. The interference between two amplitudes shown 
in Fig. ||a and b respectively generates the 4> asymme- 
try (Monte Carlo simulations with the phases of the 
bremsstralung and Ml amplitudes set so that no inter- 
ference takes place exhibit no <f> asymmetry). 

Using this model, the angular distribution in (f> is 

dT 



= Ti cos" 



T 2 sin" 



T-i sin <t> cos < 



(2) 



where the T-odd T3 sin cj) cos <f> term contains the interfer- 
ence between the Ml and bremsstrahlung amplitudes. 

Two other processes that contribute small amounts to 
the Kl tt + t:~ e + e~ decay were taken into account: the 
indirect CP-violating El photon emission (Fig. ||c) and 
the CP-conserving K° charge radius process (Fig. ||d) in 
which the Kjj — > Kg via emission of a photon. 




FIG. 3. Processes contributing to Kl — > 7r 1 7r e'e . aj 
CP-violating bremsstrahlung b) CP-conserving Ml 7 emission 
c) CP- violating El 7 emission d) Charge radius process 

The Monte Carlo simulation incorporated the ampli- 
tudes shown in Fig. ||a-d. To obtain agreement with the 
virtual photon energy spectrum E* = E e + + E e - of the 
data (Fig. [l|a), a form factor was required in the Ml vir- 
tual photon emission amplitude of Fig. ||b. We turn now 
to a detailed discussion of this form factor. 

Such a form factor has been required [p[ to explain 
the energy spectrum of the Ml photon emitted in the 
Kl — > 7r + 7r~7 decay. In order to incorporate a similar 
form factor, we have modified the coupling gyi\ 
Ml amplitude, including a form factor 

ai/a 2 



F = g M i[± 



(Ml - M\) + 2M K (E e+ +E e -Y 



(3) 



similar to that used to describe Kl — > it + tt 7 where 
M p is the mass of the p meson (770 MeV/c 2 ) and the 



photon energy has been replaced by E e + +E e ~ . The ratio 
a\jai and |<7mi| were determined by fitting the Kl — > 
7r + 7r~e + e~ data using the likelihood function 



L(ai/a 2 ,gMi) = 



H^ 1 P k M {a 1 /a 2 ,~9Mi)P^ 



(Lb Pm ( a l/ fl 2 ' 9Ml)Pa (Ol/«2, <}mi)) N 



(4) 



The probability P M of a given event is based on the 
Kl — > 7r + 7r~e + e~ matrix element and is a function of 
the five independent variables: <fi, 8 e + (the angle between 
the e + and the 7r + 7r~ 



(the angle between the ir + and the e + e 



direction in the e + e~ cms), 9 n + 
3+ °~ direction in the 
7r + 7r _ cms), M,^-, and M e + e -. It is calculated using 
the particular values of the parameters a\/a% and \§mi\ 
and nominal values from the PDG or Ref. [0 for the 
other model parameters. The likelihood of an event is the 
product of Pjjjf times P^ , the acceptance times efficiency 
of the event, normalized by the product of Pm and P a 
integrated over the entire phase space (ps). 
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FIG. 4. a) E* spectrum of data (dots), Monte Carlo using a 
constant \gMi\ (dashed histogram), Monte Carlo with E* f de- 
pendent form factor (solid histogram); b) Likelihood contours 
of ai/a,2 and |<?mi|; Constant asymmetry contours calculated 
from the data as described in the text are superimposed. 



The result of the likelihood calculation is shown in 
Fig. |b. The maximum of the likelihood occurs at 



of the a x /a 2 = -0.720 ± 0.028 GeV 2 /c 2 and \g M1 \ = 1.35 



-0.20 



where the errors represent the excursions of the likeli- 
hood function at the point where the log of the like- 
lihood has decreased by one half unit (39% CL). The 
E* spectrum predicted by the Monte Carlo with these 
parameters is shown in Fig. |]a, together with the pre- 
diction for a constant <7mi|- Figure || shows the good 
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agreement obtained using these parmeters between the 
observed <j> and sin <f> cos <p angular distributions and the 
Monte Carlo. When this form factor is included in the Ml 
amplitude, the constant |<7mi| = 0.76 used in Ref. (2| can 
no longer be directly compared to the new |<?mi| obtained 
in the likelihood fit. Rather, the average of the form fac- 
tor F of equation [3] over the range of E e + +E e - must be 
compared with the constant |<?mi| value of 0.76 ± 0.11. 
An average for F of 0.84 ± 0.10 was found, consistent 
within errors with 0.76. The branching ratio calculated 
using the form factor was increased by 5.7% compared 
with that obtained using |<?mi| = 0.76. 

Using the acceptance obtained from the Monte Carlo 
generated with the maximum likelihood values of \c/mi\ 
and ai/ci2, the asymmetry of the acceptance corrected 
sin</>cos</> distribution is found to be 13.6 ± 2.5 (stat)%. 
The contours of acceptance corrected asymmetry shown 
superimposed on the likelihood contours of a\jai and 
<7mi I in Fig. |J were determined from the sin cf> cos cf> distri- 
bution of the data, corrected for acceptances determined 
using the particular a\/a,2 and \gui\ values. 

We have considered whether the asymmetry was due to 
final state interactions. Because of the symmetry of the 
7r + 7r~e + e - state, electromagnetic or strong final state in- 
teractions, while they can modify the <j) distribution, can- 
not generate an T-odd asymmetry. The interference be- 
tween the bremsstrahlung (1=0 7r + 7r~) and the Ml (1=1 
7r + 7r~) amplitudes, which is responsible for the CP vio- 
lating asymmetry, depends on the difference of the 1=0 
and 1=1 strong interaction phase shifts but this difference 
can only modulate, not generate the asymmetry. 

Systematic errors on 01/02 and \c/mi\ due to analy- 
sis cuts, resolutions and variations of parameters of the 
Monte Carlo were studied. By varying each analysis cut 
over a reasonable range and observing the variation of 
ai/<22 and |<7mi|, ai/12 and |<7mi systematic errors of 
±0.008 GeV 2 /c 2 and ±0.04 were obtained. 

To determine the systematic errors due to resolution, 
resolution functions in the five variables were estimated 
by comparing generated and reconstructed Monte Carlo 
events. Using these functions to smear each independent 
variable for each data event, one thousand passes through 
the 1811 Kl — * 7r + 7r~e + e _ signal events were made. The 
one thousand smeared data samples were refit, and a\j a-i 
and \§mi I were determined for each of the samples. The 
variation of a\jai and \gu\\ for these samples resulted in 
errors of ±0.002 GeV 2 /c 2 and ±0.01 for ai /a 2 and \g Ml \. 

The systematic errors in a\j and \g~Mi I due to uncer- 
tainties |0] in the magnitude and phase of ry_j , and the 

uncertainties in \gEi \ and \gcR \ , estimated by varying the 
magnitude of the ratio of \gEi\ to \c/mi\ from 0.0 to 0.05 
(nominal 0.038) and \gcp\ from 0.10 to 0.17 (nominal 
0.15), resulted in systematic errors in 01/02 and gjuil of 
±0.004 GeV 2 /c 2 and ±0.01 respectively. 

All systematic errors in 01/02 and \gjsil\ were added in 
quadrature to obtain an overall error of ±0.009 GeV 2 /c 2 



and ±0.04 in 01/02 and |<7mi| respectively. 

The systematic error in the 4> asymmetry due to vari- 
ations in the corrections for acceptance arising from the 
systematic errors of the 01/02 and |<7mi and one sigma 
uncertainties of other parameters of the MC model dis- 
cussed above was determined to be ±0.7%. The variation 
in asymmetry due to analysis cuts was also estimated to 
be ±0.7%. Finally, the systematic error due to resolu- 
tion effects was determined to be ±0.7% using generated 
tracks from the Monte Carlo rather than reconstructed 
tracks in the analysis. Adding in quadrature the system- 
atic errors from these three sources, a total systematic 
error of ±1.2% was obtained for the acceptance corrected 
asymmetry of the sin (f> cos distribution. 

In conclusion, the KTeV experiment has observed a 
CP-violating asymmetry in the distribution of T-odd an- 
gle tfi in Kl — > 7r + 7r~e + e~ decays. This effect, the 
largest CP violation effect yet observed and the first 
in an angular variable, is T- violating barring possible 
exotic phenomena M such as direct CPT violation in 
the Kl — > 7r + 7r~e + e - matrix element. The magni- 
tude of the acceptance corrected asymmetry is 13.6 ± 
2.5 (stat) ± 1.2 (syst)%, consistent with the theoretically 
expected asymmetry In addition, the Ml photon 
emission amplitude requires a vector form factor as given 
in equation (3) with oi/o 2 = -0.720 ± 0.028 (stat) ± 
0.009 (syst) GeV 2 /c 2 and \g M1 \ = l-35±gf? (stat) ± 
0.04 (syst). The rich structure of the Kl — » 7r + 7r~e + e~ 
mode has provided a new opportunity for the study of 
novel CP- and T-violation effects. In the future, it may 
be possible to use this mode to search for direct CP vio- 
lation J3| and more exotic phenomena ||. 
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